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ABSTRACT 

We use N-body hydrodynamical simulations to study the structure of disks in triaxial potentials 
resembling CDM halos. Our analysis focuses on the accuracy of the dark mass distribution inferred 
from rotation curves derived from simulated long-slit spectra. We consider a massless disk embedded 
in a halo with axis ratios of 0.5: 0.6: 1.0 and with its rotation axis aligned with the minor axis of the 
halo. Closed orbits for the gaseous particles deviate from coplanar circular symmetry, resulting in a 
variety of long-slit rotation curve shapes, depending on the orientation of the disk relative to the line 
of sight. Rotation curves may thus differ significantly from the spherically- averaged circular velocity 
profile of the dark matter halo. "Solid-body" rotation curves — typically interpreted as a signature of a 
constant density core in the dark matter distribution — are obtained about 25% of the time for random 
orientations although the dark matter follows the cuspy density profile proposed by Navarro, Frenk & 
White (NFW). We conclude that the discrepancies reported between the shape of the rotation curve 
of low surface brightness galaxies and the structure of CDM halos may be resolved once the complex 
effects of halo triaxiality on the dynamics of the gas component is properly taken into account . 
Subject headings: cosmology: dark matter - galaxies: formation - galaxies: kinematics and dynamics 



1. INTRODUCTION 

It is commonly believed that the inner regions of low 
surface brightness (LSB) galaxies are ideal probes of the 
inner structure of dark matter halos. Given the small 
contribution of the baryonic component to the mass bud- 
get in these galaxies, dynamical tracers of the potential 
such as rotation curves are expected to cleanly trace the 
dark matter distribution. This provides important as- 
trophysical clues to the nature of dark matter, since the 
spatial distribution of dark material in these highly non- 
linear regions is expected to be quite sensitive to the 
physical properties of the dark matter. 

LSB rotation curves can thus be contrasted directly 
with theoretical predictions of the inner structure of ha- 
los, and there is now an extensive body of work in the 
literature that reports substantial disagreement between 
the shape of LSB rotation curves and that of circular 
velocity curves of simulated cold dark matter (CDM) 
halos (see, e.g.. iFlores fc Primackl Il994| iMoqrel Il994t 
IMcGaugh fc de Bloklll998t Ide Blok etaLl 120011) . Some 
of these rotation curves are fit better by circular veloc- 
ity curves arising from density profiles with a constant 
density "core" rather than by the "cuspy" density pro- 
files commonly used to fit the structure of CDM halos 
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(Navarro, Frenk & White 1996, 1997, hereafter NFW). 
This discrepancy adds to a growing list of concerns re- 
garding the consistency of CDM with observational con- 
straints on the scale of in dividual galaxies (see, e.g., 
iSellwood &: Kosowskvll200H) that has prompted calls for 
a radical revision of the CDM parad igm on small scales 
(see, e.g.. lSnergel fc Steinhardtil2000ft . 

Before accepting the need for radical modifications to 
CDM it is important to note a number of caveats that 
apply to the LSB rotation curve problem. For instance, 
many of the early rotation curves where the disagreement 
was noted were sig nificantly affected b y beam smear- 
ing in the HI data l|Swaters et all 12000ft . The observa- 
tional situation has now improved substantially thanks 
to higher-resolution rotation cu rves obtained from long- 
slit Ha observation s (see, e.g.. Ide Blok fc Bosmal 120021: 
LSwaters et al.ll2003ft . We shall restrict our analysis to 
these newer datasets in what follows. 

We also note that, strictly speaking, the observational 
disagreement is with the fitting formulae used to parame- 
terize the structure of simulated CDM halos (usually the 
profile proposed by NFW), rather than with the struc- 
ture of simulated halos themselves. Although the fitting 
formulae provide a simple and reasonably accurate de- 
scription of the mass profile of CDM halos, the radial 
range over which they have been validated often does 
not coincide with the scales where the disagreement has 
been identified. 

Furthermore, small but significant deviations between 
the NFW profile and simulated halos have been re- 
ported as the mass an d spatial resolution of the simu- 
lations has increased |M oore et ahl Il998t iGhigna et alJ 
l2000t iFukushige fc Makinolll997ll200H) . Although there 
is no broad consensus yet regarding how these devia- 
tions may affect the comparison with observed rotation 
curves (see, e.g.JPower et al.ll2003t lHavashi et al.ll2003t 
INavarro et alJl2004ft . the fact that the deviations worsen 
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Fig. 1. — Rotation curves of two LSB gal axies from the sa mples 
of Ide Blok fc Bosmal |2003. B02) and ISwaters et all S03), 
chosen to illustrate their v arious shapes, a s measured by the pa- 
rameter 7 from fits with the Courteau 1 1997) fitting formula. Rota- 
tion curves have been scaled to the radius ro.3 and corresponding 
velocity Vb.3 ( see text for details). Fits with large 7 values are 
characterized by a linear rise in velocity with radius followed by 
a sharp transition to the flat part of the curve. The NFW profile 
(dashed line) and the V c profiles of simulated dwarf- (dotted lines) 
and galaxy-sized halos (solid lines) match reasonably well systems 
with 7<J1 but cannot account for those with 7 > 1. 



towards the centre advise against using extrapolations of 
simple fitting formulae such as the NFW profile to assess 
consistency with observation. 

Finally, it must be emphasized that the "cusp vs. 
core" problem arises when comparing rotation speeds 
of LSB disks to spherically-averaged circular velocities 
of dark matter halos. Given that CD M halos are 
expected to be significantly non-spherical jDavis^t^d. 
1 19851 iFrenk et all 119881 Ding et all 119951 Uing fc Sutc 
2002), some differences between the two are to be ex- 
pected. It is therefore important to use the full 3D struc- 
ture of CDM halos to make predictions regarding the 
rotation curves of gaseous disks that may be compared 
directly to observation. 

We address the latter issue in this Letter, by explor- 
ing numerically the closed orbits of gaseous particles 
within the potential of an idealized triaxial CDM halo. 
We embed a massless, isothermal gaseous disk within a 
cuspy dark matter halo and evolve it until equilibrium is 
reached. 

We focus our analysis on the shape of the rotation 
curves inferred for the disk from simulated long-slit ob- 
servations of the velocity field; on deviations from the 
spherically-averaged circular velocity curve; and on the 
possibility that such deviations might account for the dis- 
crepancy between LSB rotation curves and NFW profile 
fits to CDM halos. 

2. LSB ROTATION CURVES 

Figure illustrates the disagreement alluded to above. 
This figure shows the Ha rotation curves of two LSB 
galaxies (points with error bars) selected from the sam- 



ple of Ide Blok fc Bosmal <|2002l B02) and iSwaters et alJ 
(2003, S03). The data points have been fitted using a 
simpl e formula, V TO t(r) = Vo(l+(r/r t ) -7 ) -1 / 7 ijCourteaul 
I1997J) . Here Vq and r t are dimensional scaling parame- 
ters, whereas 7 is a dimensionless parameter that char- 
acterizes the shape of the rotation curve. This three- 
parameter formula provides excellent fits to both LSBs, 
as illustrated by the quality of the (solid line) fits shown 
in Figure ^ 

In order to emphasize discrepancies in shape, the ro- 
tation curves in FigureQhave been scaled to the radius, 
ro.3, and velocity, V0.3, where the logarithmic slope of 
the curve is d log V ro t /d log r = 0.3. The two galaxies 
shown in Figure ^ have different values of 7, and have 
been chosen to illustrate the extreme cases in the B02 
and S03 datasets. Roughly one third of their LSBs have 
7<1, having rotation curves with shapes similar to UGC 
5721 (7 = 0.6). 

The dashed line in Figure ^ shows the V c profile of an 
NFW halo, which is fixed in these scaled units. Galax- 
ies with 7<1 are consistent with NFW, whereas those 
with 7 ^> 1 are clearly inconsistent. Figure H also shows 
the spherically- averaged V r profiles of all g alaxy-sized ha- 
los pre sented in lHavashi et al.l l)2003j) and lNavarro et all 
(2004), scaled to ro.3 and V0.3. We find that the shapes 
of the dark halo V c curves are in fact quite similar to 
NFW. In terms of the 7 parameter, most halos (about 
95%) have 7<1, which implies that LSBs with 7^1 are 
quite difficult to reconcile with the V c profiles of simu- 
lated CDM halos. 

On further examination, however, lHavashi et alJ 
( 2003) note that most rotation curves that have best fit 
values of 7 > 1 also have acceptable fits with 7 < 1. As 
a result, only a small minority of LSBs (about 10%) are 
robustly inconsistent with CDM halo V c profiles. Most of 
these curves are characterized by a linear rise in velocity 
with radius and have best fit values of 7>5. One such 
example, UGC 5750, is shown in Figure ^ along with 
a dot-dashed line that illustrates the V oc r dependence 
expected in the presence of a constant density core. 

Does this result rule out the presence of a cusp in the 
dark matter density profile in such galaxies? As noted in 
§n before concluding so one must take into account pos- 
sible systematic differences between rotation speed and 
circular velocity in gaseous disks embedded within real- 
istic, triaxial halos. This is a complex issue that involves 
a number of parameters, such as the degree of triaxiality, 
the role of the disk's self-gravity, size, and orientation, as 
well as the possibility of transient deviations from equi- 
librium. 

To keep matters simple, we have decided to ad- 
dress this issue by evolving a massless gaseous disk 
at the centre of a fixed triaxial halo. We use the 
N-body/hydrodynamical code GA SOLINE, developed b y 
J.Wadsley, J.Stadel, and T.Quinn l)Wadslev et all 2004). 
GASOLINE combines a tree-based Poisson solver for gravi- 
tational interactions with the Smooth Particle Hydrody- 
namics (SPH) technique. The dark matter halo is mod- 
elled with a particle realization of an NFW mass profile; 
the mass, virial radius and concentration of the halo are 
normalized to Af2oo = I0 11 M & /h, r2oo = 75 /i _1 kpc and 
c = 12, respectively. In order to reduce computational 
expense, the halo is truncated at an outermost radius, 
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Fig. 2. — Evolution of a massless gaseous disk in the 3D poten- 
tial of an idealized CDM halo. Upper (lower) panel shows edge-on 
(face-on) projections of the disk. Left panels shows halo and disk 
particles at the initial time, color-coded by local dark matter and 
gas density, respectively. Middle and right panels show disk parti- 
cles only at initial time and after 7.0 Gyr, respectively. Inner and 
outer circles have radii equal to the exponential scale length of the 
disk, r,j, and the outermost radius of the disk, r ou tcr, at t = 0, re- 
spectively. The straight lines in the right panels show projections 
of the line-of-sight used to generate the rotation curve in the lower 
right panel of Figure [3] Triaxiality in the halo mass distribution 
leads to significant evolution in the structure of the disk and to 
strong deviations from circular symmetry. 



''trunc = 37.5 /i _1 kpc and the number of particles within 
this radius is set to 1.7 x 10 5 . The halo is made triaxial 
by multiplying the x-, y-, and z-coordinates of the halo 
particles by factors of 1.36, 0.82, and 0.68, respectively, 
so that the spherically-averaged mass profile remains un- 
changed, but the axis ratios of the halo are 0.5: 0.6: 1.0. 
We note that such halo shapes (roughly prolate with 
elongation 2:1) are among the most common ones foun d 
in cosmological simulations (see, e.g.. Uing fc Sutol l2002). 
The gravitational potential is derived directly from the 
dark matter particle positions and is assumed to remain 
constant in time. 

An exponential disk consisting of 10 4 massless gas par- 
ticles is placed at the centre of the halo, in the plane 
of the major and intermediate axes of the halo (see 
Figure |2J) ■ The exponential scale length of the disk is 
r d = 3.6 /i _1 kpc, and its outermost radius IS Pouter — 
18.0 /i _1 kpc. The vertical distribution of the disk mate- 
rial is given by p(z) cx sech 2 (z/zo), with zq — 0.6 /i _1 kpc, 
and the disk is truncated at z max = 3.0 ft kpc. The 
gas is given a temperature T = 100 K, well below the 
virial temperature of the halo, T V j r ~ 2 x 10 5 K, and 
is assumed to remain isothermal during the evolution. 
Pressure forces are thus unimportant for the dynamics 
of the gas, and the hydrodynamical treatment simply 
forces the gas to provide a massless tracer of the closed 
orbits within the halo potential. 

The disk particles are initially given tangential veloci- 
ties consistent with the spherically-averaged circular ve- 
locity profile of the triaxial dark matter halo. As a re- 
sult the gaseous disk is not initially in equilibrium and it 
evolves rapidly to a configuration characterized by depar- 
tures from circular symmetry (the disk becomes ellipti- 
cal, see Figure[3J|. The effect of such departures from cir- 
cular symmetry on rotation curves derived from long-slit 
spectra is complex, as illustrated in Figure^ This figure 
compares the rotation speed of the g ciS, clS inferred from 
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Fig. 3. — Upper panels: Projected positions of the disk particles 
after the disk has been inclined by 67° relative to the initial plane 
of the disk. Disk particles are color-coded by line-of-sight velocity. 
Solid horizontal lines indicate the position of a 2 h~ 1 kpc wide 
slit oriented along the photometric major axis of the projected 
disk. Lower panels: Rotation curve as inferred from simulated 
long-slit radial velocity data with slit placed across the major axis 
(open circles) and minor axis (open triangles) of the disk. The 
major axis rotation curve agrees well initially (t = 0) with the V c 
profile of the halo (dotted curves) but significant deviations and 
misalignment of the kinematic and photometric major axes result 
from the evolution of the disk in the triaxial potential of the halo. 



line-of-sight velocities measured on a slit placed along 
the photometric major axis of the disk determined from 
the projected particle distribution. The disk is inclined 
by 67°, and velocities are corrected by the sine of the 
inclination angle, z b s , as derived from the aspect ratio 
of the isodensity contours of the gas. 1 

The left panels in Figure|21show that initially, when the 
disk is circularly symmetric by construction, the disk in- 
clination is recovered reasonably accurately and the ro- 
tation speed inferred from line-of-sight velocities (open 
circles) agree well with the spherically-averaged circular 
velocity (dashed line). Fitting the halo V c profile out to 
foutcr = 18 /i kpc using the (r t , 7, Vq) formula results in 
best fit values of 7 = 1, whereas the initial disk rotation 
curve is best fit by 7 = 1.3. 

At later times the evolution of the disk leads to poorer 
estimates of the inclination as well as to significant de- 
viations between inferred rotation speed and circular ve- 
locity. The shape of the rotation curve, in particular, is 
affected, as shown in the right-hand panels in Figure [3J 
On some projections, rotation curves appear to rise and 
turn abruptly, and they would be (erroneously) taken to 
imply the existence of a constant-density core in simple 
models that assume spherical symmetry. Fits to the disk 
rotation curves using the (rt,7, Vo) formula introduced 
above often have 7 > 1, consistent with galaxies where 
NFW profiles provide a particularly poor fit to the rota- 
tion curve data. The open triangles in the lower panels 

This is inferred to be ~ 58° because the disk is not circularly 
symmetric. 
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Fig. 4. — Systematic deviations from circular velocity induced by 
the triaxiality of the halo lead to rotation curve shapes in better 
agreement with LSBs with high 7. Open circles correspond to 
rotation curves inferred from simulated long-slit data for the disk 
shown in Figure |3 



of Figure |3| show the rotation curves obtained with the 
slit oriented along the minor axis of the disk. The mis- 
alignment between the photometric and kinematic axes 
of the disk after 7.0 Gyr is clearly evident. 

We have calculated the distribution of 7 values for ro- 
tation curves obtained by observing the simulated disk at 
time t = 7.0 Gyr from 1000 random lines-of-sight with 
inclination angles limited to 30° < i < 70°. A large 
fraction of rotation curves have values of 7 significantly 
higher than the initial value of 7 = 1.3. To be precise, 
approximately 45%, 25%, and 7%, have 7 > 2, > 3, and 
> 4, respectively. Furthermore, minor axis rotation is 
seen with a peak amplitude exceeding 20% of the major 
axis amplitude in 70% of simulated rotation curves with 
7 > 3. 

In Figure^! we compare the high-7 LSB rotation curve 
shown in Figure^with one obtained from the projection 



of the disk shown in Figure The agreement between 
the simulated and observed rotation curves is excellent, 
suggesting that deviations from spherical symmetry in 
CDM halos might reconcile disk rotation curves that ap- 
pear to favour the presence of constant density cores with 
cusps in the dark matter density profiles. 



3. DISCUSSION 

Given that we are able to match discrepant rotation 
curve shapes with a disk embedded in a cuspy triax- 
ial halo, the outlook for reconciling dark matter cusps 
with LSB rotation curves is rather encouraging. How- 
ever, it would be premature to argue that the problem 
has been fully solved. After all, given the number of extra 
"free" parameters introduced by relaxing the assumption 
of spherical symmetry, it is perhaps not surprising that 
one is able to improve the agreement with LSB rotation 
curves. 

It is therefore important to build a more compelling 
case for this interpretation of LSB rotation curves, so as 
to render it falsifiable. Are there any corroborating traits 
that may be used to confirm or exclude the hypothesis 
that halos surrounding LSBs are indeed triaxial? In par- 
ticular, we would like to understand better the particular 
combination of perspective and triaxiality that results in 
rotation curves with values of 7 > 1. How can one best 
verify the triaxial-halo interpretation in two-dimensional 
velocity maps? Identifying a clean and unambiguous in- 
dication of triaxiality, such as the unusual minor axis 
kinematics shown in Figure will be as important as 
the success of aspherical halos in reproducing the rich 
variety of shapes of LSB rotation curves. Only if this is 
accomplished shall we be able to conclude that LSB rota- 
tion curves do not preclude the presence of dark matter 
density cusps, thereby freeing the CDM paradigm of one 
vexing challenge on small scales. 
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